Scattering of surface plasmons by one-dimensional periodic nanoindented surfaces 
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In this work, the scattering of surface plasmons by a finite periodic array of one-dimensional 
grooves is theoretically analyzed by means of a modal expansion technique. We have found that 
the geometrical parameters of the array can be properly tuned to achieve optimal performance of 
the structure either as a Bragg reflector or as a converter of surface plasmons into light. In this last 
case, the emitted light is collimated within a few degrees cone. Importantly, we also show that a 
small number of indentations in the array are sufficient to fully achieve its functional capabilities. 
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Surface plasmons (SPs) are well known for their ca- 
pabilities for concentrating light in sub-wavelength vol- 
umes and guiding light through the surface of a metal 
1]. This has recently raised the prospect of SP-based 
photonic circuits, with length scales much smaller than 
those currently achieved 0,0,0. In order to reach this 
goal, it is mandatory to know the scattering properties of 
SPs by simple surfaceprofiles, aspointed out by recent 
experimental studies [3, la lj, Ld- lj- III From the theoretical 
side, this is still an open problem, even after a few sem- 
inal works based on the Rayleigh expansion [lH ITU ll2T ] 
and the Green's function dyadic 0, Q|. Both these 
methods require the evaluation of sophisticated scatter- 
ing functions from which physical insight is not easily 
inferred. 

In this work we present an alternative formalism for 
calculating the scattering of SPs by indentations per- 
forated on a thick metal film. For that purpose, we 
have extended to real metals the modal expansion tech- 
nique previously developed within the perfect conductor 
approximation (PC A) [id llf| . therefore incorporating 
surface plasmon polaritons into the model. This ap- 
proach enables us to describe the scattering properties 
of an arbitrary set of indentations without any restric- 
tion over their position or shape. Besides, our method 
does not require any adjustable parameter, inasmuch as 
the wavelength-dependent dielectric function of the metal 
e(\) is previously known. In addition to this, it provides a 
very compact representation of the electromagnetic (EM) 
fields and simple expressions for the scattering magni- 
tudes, which permits us to extract underlying physical 
mechanisms much more easily. 

Within our theoretical framework, the way of launch- 
ing SPs onto a set of indentations resembles the back-side 
illumination employed in some experimental works 
We consider a single slit flanked by a set of N indenta- 
tions placed in the output surface of an infinite metallic 
film of thickness h (see Fig.l). Eventually, the distance 
between the slit and indentations will be taken to be 
infinity. In this way, the slit merely plays the role of 



a theorist's SP-launcher. More precisely speaking, SPs 
are taken into account by applying surface impedance 
boundary conditions (SIBC) 17] to the metal/dielectric 
interface along the film surface. 

In this work, we consider the simplest case of ID sub- 
wavelength indentations (grooves). Additionally, we im- 
pose that the external illumination be uniform along the 
y-axis, so we restrict ourselves to scattering of SPs im- 
pinging onto the grooves at normal incidence, where only 
the fundamental mode inside each groove is relevant. 
This leads to a set of N + 2 equations for the modal 
amplitudes of the electric field at the input and output 
openings of the slit (E , E' Q ) and the output openings of 
the N grooves ({££}): 

(Goo — eoj-E-o — GvoE = la- 
(G aa — e a )E' a + ^^GapE'p — GyoEoSaQ = 0. 

(1) 

Notice that this set of equations is the same as the one 
obtained within the PC A Therefore, its physi- 

cal meaning remains unchanged, although the expres- 




FIG. 1: Schematic picture of the system under study: single 
slit flanked in the output surface by an arbitrary set of N 
indentations located at a long distance from its right side. A 
p-polarized EM wave is impinging from the bottom. Param- 
eters {x a ,a a ,w a } define the geometry of indentations. 
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sions are slightly different due to the non-zero impedance 
Z s = e(A) -1 / 2 at the metal surface: Iq takes into account 
the back-side illumination on the slit (which must be p- 
polarized in order to excite SPs); the "self-energy" e Q is 
related to the bouncing back and forth of the EM fields 
inside indentation a. For a groove, e a = —i(l + (f> a )/(l — 
Z s -(l + Z s )c(> a ), where 4> a = e 2ft « w "(l - Z,)/(l + Z s ), 
being ko = 2tt/X and w a the depth of the groove; 
G v0 = 2i/(e zkoh (l + Z s ) 2 - e - lkah {l - Z s ) 2 ) represents 
the coupling between the two sides of the slit. Finally, 
Ga/3 is the coupling between modes, reflecting that mode 
(3 emits radiation that can be collected by mode a. More 
precisely, G a p is the projection onto the wavefields at 
the openings of indentations a and f3 of a scalar Green's 
function 



G{x^ x , z ) 



A 



dk- 



i(k\x-x'\ + y/k%-k 2 \z-z'\) 



yjkl -k 2 + k Z s 



(2) 



As Z s — > 0, Eq. (j2J) transforms into an integral represen- 
tation of the Oth-order Hankel function of the first kind 
(lif. thus recovering the PCA in the low-energy limit of 
our approach. Even within the SIBC approximation, we 



find that the PCA result is still valid for 



<< A. 



However, the presence of Z s in Eq.J2J strongly modifies 
its long-distance behavior. By means of the saddle-point 
approximation, we have obtained an asymptotic expan- 
sion of G in the limit where \z — z'\, A <C \x — x'\: 

G as (x, z; x 1 , z') = -(k 2 Z s /k p ) e i(fc,l*-*'l-fcaZ.I*-*'l) ) 

(3) 

with k p satisfying the SP dispersion relation of a flat 
metal-dielectric interface within the SIBC 



kp — 



— Z s kn. Therefore, the long distance EM coupling along 
the surface is due to SPs, even in the presence of ab- 
sorption. Comparison with the exact result in the op- 
tical regime shows that the asymptotic limit is already 
reached for small distances. For example, in the case of 
silver at A = 750 nm, we find that G as (f, r') differs from 
G{f,r') by less than 10% for \x - x'\ w 2A. It is this 
knowledge of long-distance EM coupling being mediated 
by plasmons what allows us to use the system in Fig. ^ 
for the analysis of SP scattering. 

Let us then examine the term Go a E' in Eq. Q. It 
can be interpreted as an "slit illumination" impinging on 
the grooves. Thus, the equation governing the EM fields 
at the grooves becomes 



(G a a 



Gj3 a E'p — I a 



(4) 



where I a = —Go a E' is defined to resemble the back-side 
illumination Iq. The key point is that, according to Eq. 
0, I a corresponds to a SP illumination, modulated by 
a constant factor that depends on the metal thickness, 
the intensity of back illumination and the width of the 



slit. This factor is not relevant for the determination 
of scattering coefficients and the whole slit can then be 
treated as a theoretical artifact. 

Once the self-consistent {E' a } are obtained, the calcu- 
lation of the EM field in all space is straightforward, and 
therefore both the emittance S (which is the fraction of 
incoming SP energy radiated into vacuum) and its angu- 
lar distribution. As the EM coupling between the grooves 
and a distant point on the surface is due to SP, we can 
also obtain the SP reflection (r) and transmission (t) am- 
plitudes: 
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c a e ikpXc 

a=l 
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(5) 



where c a = — ifcg V / a^(Z s /fc p )sinc (k p a a /2) is a geomet- 
rical coefficient associated to each indentation. Notice 
that, if absorption is present, im[fc p ] ^ and the SP 
reflected and transmitted currents depend on the points 
(f R , r T ) at which the EM are evaluated, reflecting the ab- 
sorption loss in the flat regions of the metal surface. This 
suggests that the scattering coefficients should be ex- 
tracted from the EM fields at points close to the grooves, 
although in this case it is difficult to separate the diffrac- 
tive contribution from the one due to SPs. Nevertheless, 
provided that the grating lengths considered are shorter 
than the SP absorption length, absorption can be ne- 
glected. In what follows, we present the results obtained 
under such assumption for finite periodic arrays of N 
rectangular grooves, patterned on a Ag film. 

We consider grooves with width a — 100 nm and depth 
w = 50 nm separated by a period d of 600 nm, which are 
typical experimental values. Figures H[a)-(c) render the 
calculated reflectance R = |r| 2 , transmittance T = \t\ 2 
and emittance S spectra, for increasing values of N. For 
a single groove (top panel), T increases with A, while both 
R and S decrease. This is due to two mechanisms. First, 
there is a decrease of the relative size of the groove with 
respect to A, which manifests in G a p scaling as (a/ A). 
Second, at longer As, the SP wavefield is more extended in 
the air region and therefore less sensitive to the presence 
of obstacles at the surface. Panels (b) and (c) show how 
the addition of more grooves greatly modifies the optical 
response of the system. As N increases, transmission 
gaps develop, as well as sharp resonances in both R and 
S. 

In order to gain insight into the origin of this behav- 
ior, it is helpful to analyze the EM surface modes of an 
infinite groove array. This can be readily done by look- 
ing for solutions to Eq. @, imposing both I a = and 
Bloch's theorem (i.e. E' a = E'e lk * ad , k x being the ex- 
tended surface state wavevector at the given wavelength). 
The band structure (solid line) for surface modes in a pe- 
riodic structure with the same geometrical parameters as 
in (a) to (c) is presented in Fig^d), as well as the disper- 
sion relation of SPs in a flat air/silver interface (dots). As 
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expected |l9|. band gaps occur with a low-A edge given 
by k p d = mm with m = 1, 2, • • • , i.e. by the folding of the 
dispersion relation in a flat surface. On the contrary, the 
high-A edge depends on the geometry of the grooves, as 
it corresponds to a SP standing wave with maxima at the 
groove positions. Evidently, spectral regions of low T in 
the finite array coincide with gaps in the band structure. 
Energy conservation implies a corresponding increase in 
R + S, but it is not obvious how this increase is divided 
between these two channels. However, there is a simple 
argument for the existence of reflection maxima. Let us 
consider the SP wavefields emitted by two consecutive 
grooves in the region of reflection. There is an "optical 
path" phase difference of k p d between these waves. Ad- 
ditionally, there is also a phase difference between emit- 
ters that is equal to k x d in the case of a infinite system. 
But, as previously noted, k x d — k p d — mir at the low- 
wavelength gap edge so the SP wavefields launched by all 
grooves interfere constructively (notice that they also in- 
terfere constructively in the transmission region, but not 
with the incident field) . As A is increased away from this 
condition, the constructive interference is progressively 
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FIG. 2: (color online). Calculated R (solid), T (dashed) and 
S (dotted) curves of SP scattering by a finite periodic groove 
array on a Ag film. Here, a = 100 nm, d = 600 nm and w = 50 
nm. Results in (a), (b) and (c) correspond to structures with 
1, 5 and 10 grooves respectively. Panel (d) shows the band 
structure (solid lines) for the same parameters as in (a) to 
(c) and the SP dispersion relation in a flat air/Ag interface 
(dots). Gray strips mark the photonic gaps. 
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FIG. 3: (color online). Maximum values of R (a) and S (b) as 
a function of the number of grooves in the vicinity of k p d = 2ir 
for the same w, d as in Fig. Hand increasing values of a. 



lost, R decreases and, for A within the gap, S increases. 
If A crosses the gap edge, the transmission channel is 
open, and S decreases. Therefore, S presents peaks at 
the high-A edges of the gaps, as can be seen in Fig|5] 
In our opinion, this mechanism can explain the heuristic 
criterion for optimum mirror efficiency presented in Ref. 
and may also be at the root of the strong asymmetry 
in the positions of the reflectance peaks reported in Ref. 
i 

Figure renders the maximum R and S as a function 
of the number of grooves in the vicinity of k p d = 2tt 
for the same w,d parameters as in Fig. [2 and increasing 
values of a. Notably, a small number of indentations are 
sufficient to achieve either a large in-plane reflection or a 
high emission out of the plane. This rapid saturation is 
also consistent with Ref. 6] . With respect to the groove 
width, it is clear that it mainly influences the out-of-plane 
efficiency of every single scatterer, being more relevant for 
the S vs T ratio than for mirror efficiency. 

The dependence of the scattering coefficients on groove 
depth is depicted in Figure for N — 10, a — 100 nm 
and d = 600 nm. As can be seen, the maxima of R vary 
weakly on w for most ranges of groove depth. Such a 
weak dependence may be relevant for device design, con- 
sidering that the control of w is often the most difficult 
point in groove fabrication. However, Figure0|also shows 
that for some values of w the reflection is very small, 
when a maximum is expected. Simultaneously, at these 
groove depths, S sw and T w 1. This occurs close to 
the condition A = 4w/ (2n + l), when the in-plane electric 
field at the indentations is very small, due to destructive 
interference of the incident field and the field reflected on 
the closed end of the groove. 

One of the possible applications of finite arrays of in- 
dentations lies in their capability to convert SPs into 
light. Therefore, it is worth studying the directional- 
ity properties of the emitted light in the system analyzed 
throughout this paper. Figure shows the far-field an- 
gular distribution of radiation emitted out of plane, eval- 
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FIG. 4: (color). Contour plots of R (a), T (b) and S (c) 
versus both groove depth and wavelength for the same N, a, d 
as in Fig. Etc). 

uated for the two emittance maxima in Fig. 2c (A = 665 
nm and A = 500 nm) at the low-energy edges of the gaps 
labelled, respectively, with m = 2 and m = 3. The dis- 
tribution corresponding to the gap with m = 2 (solid 
curve) is beamed close to the normal. On the contrary, 
the one coming from the gap with m = 3 is beamed at 
higher angles. Notice that, at the condition of maximum 
R, E a cx (— l)( ma ), so the emittance would be normal for 
even m and close to tangent for odd m. Inboxes render 
the FWHM of angular distribution peaks for both curves 
as a function of the number of grooves. As expected from 
usual grating theory, it scales as 1/N . 

To summarize, in this paper we have presented a theo- 
retical framework that can treat the scattering of SPs 
by a finite array of indentations. With this formal- 
ism, we have studied the scattering properties of a sub- 
wavelength periodic groove array in the visible and near 



IR ranges. We have found that, associated to the low-A 
edge of the SP band gap, the array behaves as a mir- 
ror (up to 80% reflectance) whereas at the high-A edge, 
most of the light carried out by the SP can be converted 
into collimated light (up to 90%). We have also shown 
that this resonant behavior can be achieved with a small 
number of indentations (around 10) and it is quite robust 
with respect to variations in the groove depth. 
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